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ABSTRACT 



Context. Open clusters offer a unique possibility to study the time evolution of the radial metallicity gradients of several elements in 
our Galaxy, because they span large intervals in age and Galactocentric distance, and both quantities can be more accurately derived 
than for field stars. 

Aims. We re-address the issue of the Galactic metallicity gradient and its time evolution by comparing the empirical gradients traced 
by a sample of 45 open clusters with a chemical evolution model of the Galaxy. 

Methods. At variance with previous similar studies, we have collected from the literature only abundances derived from high- 
resolution spectra. The clusters have Galactocentric distances 7 < R GC < 22 kpc and ages from ~ 30 Myr to 11 Gyr. We also 
consider the ff-elements Si, Ca, Ti, and the iron-peak elements Cr and Ni. Cepheids trace instead the present-day Fe gradient in the 
inner parts of the disk. 

Results. The data for iron-peak and or-elements indicate a steep metallicity gradient for R GC < 12 kpc and a plateau at larger radii. 
The time evolution of the metallicity distribution is characterized by a uniform increase of the metallicity at all radii, preserving 
the shape of the gradient, with marginal evidence for a flattening of the gradient with time in the radial range 7-12 kpc. Our model 
is able to reproduce the main features of the metallicity gradient and its evolution with an infall law exponentially decreasing with 
radius and with a collapse time scale of the order of 8 Gyr at the solar radius. This results in a rapid collapse in the inner regions, i.e. 
^gc ;S 12 kpc (that we associate with an early phase of disk formation from the collapse of the halo) and in a slow inflow of material 
per unit area in the outer regions at a constant rate with time (that we associate with accretion from the intergalactic medium). An 
additional uniform inflow per unit disk area would help to better reproduce the metallicity plateau at large Galactocentric radii, but it 
is difficult to reconcile with the present-day radial behaviour of the star formation rate. 

Conclusions. Our results favour a scenario where the Galactic disk is formed inside-out by the rapid collapse of the halo and by a 
subsequent continuous accretion of intergalactic gas 

Key words. Galaxy:disk, abundances, evolution; open clusters and associations: general, abundances 



1. Introduction 

The radial metallicity gradient and its time evolution provide 
strong constraints on our understanding of the formation and 
the evolution of galaxies. Detailed models of Galactic chemical 
evolution (GCE) have been continuously improved over the last 
decad es (e.g., Lacey & F all[l985l Tosi [l988l Ferrini et al. [l992l 
I1994I Giovagnoli & Tosi |1995l Chiappini, Matteucci & Gratton 
[19971 Boiss ier & Prantzos [1999t Hou et al. |2000l Chiappini 
et al. 12.001b . However, the issue of the time evolution of the 
radial abundance gradients is far from being settled unequiv- 
ocally, both from a theoretical (e.g., Goetz & Koeppen]1992; 
Koeppen [19941 Molla et al. [19971 Henry & Worthey 1999; 
Chiappini et al. 120011 ) and an observational (e.g., Friel 1995; 
Friel et al. 120021 Maciel |2000l Maciel et al. |2U03l Stanghellini 
et al. 120061) point of view. While most GCE models are able to 
reproduce the present-day radial distribution of several chemi- 
cal elements derived from sample objects representative of the 
present-day composition of the interstellar medium (ISM), such 
as, e.g., H n regions, Cepheids, super-giant stars, they generally 
disagree on the predicted behaviour of its time evolution. 

In this respect, GCE models can be divided in two groups: 
models which predict a steepening of the metallicity gradient 
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wi th time (e.g., Tosi 119881 Chiappini et al. |1997t Chiappini et 
al. 2001) and those which predict a flattening with time (e.g., 
Molla et al. [19971 Portinari & Chiosi [19991 Hou et al.[2000). 
The main differences between the two groups of models are the 
efficiency of the enrichment processes in the inner andf outer 
regions of the Galactic disk (the star formation rate, SFR), and 
the nature of the material (primordial or pre-enriched) falling 
from the halo onto the disk (the infall). 

In models of the former group the outer disk is pre-enriched 
by the previous evolution of the halo, and during the first Gyr its 
metallicity is affected very little by the relatively lower SFR in 
the disk. Since the infall decreases rapidly with Galactocentric 
radius, reflecting the strongly concentrated mass distribution in 
the halo, a larger amount of metal poor material falls in the cen- 
tral regions rather than in the outer Galactic disk. Thus, a pos- 
itive metallicity gradient is initially established. When the halo 
collapse phase is over, the metallicity increases more rapidly in 
the inner disk than in the outer Galaxy, because of a combination 
of higher SFR at small radii and continuous dilution by infall at 
large radii, reversing the sign of the gradient. 

In contrast, in models of the latter group, the initially vig- 
orous star formation activity in the inner disk and the metal- 
poor composition of the infalling gas produce a rapid increase of 
the metallicity near the Galactic center. In the inner regions the 



2 



Magrini et al.: Evolution of the Galactic metallicity gradient 



metallicity increases rapidly, reaching its final value in the first 
2-3 Gyr of disk evolution, whereas the enrichment of the outer 
disk is slower, resulting in a progressive flattening of the radial 
gradient. In this group of models, various infall rates have been 
assumed: a very rapid one simulating the formation of the disk 
from a monolithic collapse of the halo, or one much diluted with 
time, simulating a hierarchical formation of the disk by contin- 
uous infall of gas from the intergalactic medium (IGM). These 
rates produce different time evolution of the metallicity gradi- 
ent: a rapid flattening or a uniform increase of the metallicity 
with essentially no change of slope, respectively (e.g., Molla et 
al JT996l 

In principle, it should be possible to discriminate between 
the different evolutionary scenario, by comparing the metal- 
licity gradient of the "old" stellar population e.g., red giant 
branch (RGB) stars, planetary nebulae (PNe), and old open clus- 
ters (OCs), with the present-day gradient outlined, e.g., by H n 
regions, young stars (including Cepheids and giant stars) and 
young OCs. In external galaxies, the comparison of PN and H n 
region abundances, investigated with similar observational and 
analysis techniques, has given encouraging results (e.g., Magrini 
et al. |2007| for the spiral galaxy M33), since PNe can be assumed 
at the same distance as the host galaxy. 

In our Galaxy, however, the huge uncertainty on the dis- 
tances to Galactic PNe restrict their use as reliable tracers of 
the past evolution of metallicity gradients. On the contrary, the 
determination of metallicities of OCs represents a reliable ap- 
proach to the study of the time evolution of the metallicity gra- 
dient. OCs represent, in fact, a unique opportunity to study abun- 
dances of a large number of elements in stars of known ages and 
located at a wide range of well-determined Galactocentric dis- 
tances. Thus, they can fully address the issue of the time evo- 
lution of the metallicity gradient, provided that a homogeneous 
analysis is performed. 

Several spectroscopic studies of chemical abundances in 
OCs have been carried out during the last years, both at low- 
(e.g., Friel [19951 Fri el et al. |20U2l and high-re s olutio n (e.g., 
Carretta et al. [20041 Carrar o et al. 12004, 2006J |2007l Yong, 
Carn ey & d e Almeida 120531 Sestito et al l2006ll2008l Bragaglia 
et al. 120081 ). Different authors find different results, depending 
on the quality of the data, and in particular on the spectral reso- 
lution, and on the method of analysis adopted. These discrepan- 
cies especially affect metallicity gradients. Also, abundances for 
the farthest clusters, which are critical to determine the gradi- 
ent up to very large Galactocentric distances, were not available 
until very recently. For example, Friel ( 119951 [2006), Carraro, 
Ng & Portinari ( 119981 ) and Friel et al. ( 120021 ) suggest the pres- 
ence of a negative [Fe/H] gradient, while Twarog, Ashman & 
Anthony-Twarog (1997) and Corder & Twarog (200T) favor a 
step-like distribution of the Fe content with Galactocentric dis- 
tance. Furthermore, recent results (e.g., Yong et al. 120051 Carraro 
et al. |20071 Sestito et al.|2007) suggest that the gradient becomes 
flat for Galactocentric radii Rcc larger than ~ 11 — 12 kpc. 

For the time evolution, evidence for variation so far has been 
not conclusive. The majority of studies seem to favor a slight 
evolution from steeper gradients in the past to shallower one at 
the present time. For example, Friel et al. (120021) . based on metal- 
licities derived from low resolution spectra, find a flattening of 
the overall gradient from -0.096 dexkpc 1 for ages greater than 
6 Gyr, to -0.023 dexkpc 1 for clusters younger than 2 Gyr; sim- 
ilarly, Chen et al. ( 120031 ) derive a shallower slope for clusters 
younger than 0.8 Gyr than for older ones. On the contrary Salaris 
et al. ( 120041 ) find a steeper slope for younger clusters. These dis- 
crepant results are likely due to: i) use of different and/or inho- 



mogeneous abundance datasets; ii) the different radial regions of 
our Galaxy which are analyzed. 

In this paper we re-address the issue of the radial metal- 
licity gradient in our Galaxy observationally and theoretically, 
exploiting the availability of abundances based on high reso- 
lution spectra for several open clusters. Specifically, we col- 
lected and analyzed the most recent metallicity determinations 
of a large sample of OCs, considering only measurements based 
on high resolution spectroscopy, which by far are more reliable 
than those based on photometry or low resolution spectroscopy. 
These observations, described in Sect. [2l are compared with the 
results of our chemical evolution model, able to reproduce the 
main observational features of our Galaxy. In addition to Fe, the 
most recent studies (e.g., Friel et al. 120031 Carretta et al. 2004; 
Friel Jaco bson, & Pilacho wski 120531 Yong et al. 120051 Sestito 
et al. 2006J Randi ch et al. [20061 Bragaglia et al. [20081 Sestito 
et al. 120071 120081 ) also include the analysis of other elements. 
These studies allow us to investigate for the first time the radial 
gradient of Fe-peak and a-elements based on the comparison of 
OC abundances, rather than on field stars, and model predictions. 
Finally, we consider the empirical gradient traced by Cepheids 
for the inner parts of the disk (Rgc ^ 7 kpc). 

The rest of the paper is organized as follows: in Sect. [3] we 
describe the model; in Sect. [4] we analyze the time evolution of 
the radial gradient of several chemical elements, including Fe 
and several a and iron-peak elements; the abundance ratio radial 
behaviours are discussed in Sect. 14.21 while the results of this 
study are discussed in Sect.[6l in Sect. [7] we summarize our con- 
clusions. In the Appendix we describe the details of the Galactic 
chemical evolution model. 

2. The data sample 

As discussed in Section 1, since the early study by Janes ( 1979) 
several works have focused on the determination of the radial 
metallicity gradient based on OCs. The first remarkable studies 
addressing the time evolution of the gradient through OC data 
were performed by Friel et al. (2002) and by Chen et al. (2003 ). 
The former analyzed a sample of low- and moderate-resolution 
data of Galactic OCs over a range in Galactocentric radius from 
7 to 16 kpc. The latter was based on a collection of metallicities 
for a very large sample of OCs up to a distance <16 kpc from 
the Galactic center, whose metallicities and ages had been de- 
rived by different authors and with different methods. Therefore, 
whereas the catalog of Chen et al. (2003) represents the largest 
sample of OCs so far collected to study the gradient, it is greatly 
inhomogeneous. 

In the last few years, and after the two studies mentioned 
above, thanks to state-of-the-art multiplex facilities on 8m class, 
the number of OC with secure metallicity determinations from 
high-resolution spectroscopy has greatly increased. Also, and 
most important, it has been possible to extend metallicity deter- 
minations to OCs in the outer disk, beyond 15-16 kpc from the 
Galactic center. Finally, abundances of other elements besides 
iron have been determined. 

In this context, the strength of our data sample is that it com- 
bines the largest available sample of high-resolution observa- 
tions of OCs, including OCs with Galactocentric distances be- 
tween ~ 7 and 22 kpc, together with a homogeneous treatment 
of their age derivation, attaining larger Galactocentric distances 
than in previous samples. 

Our dataset of OCs (see Table Q]) is a collection of sam- 
ples analyzed by different authors in the literature. Table Q] lists 
in Cols. 1 and 2 the sample clusters and their Galactocentric 
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distances Rgc, taken from Friel et al. (|2002| l, when available, 
otherwise from Friel ( 119951 120061 ); in the few cases when Rgc 
from Friel was not available, we calculated this quantity using 
the cluster distance given in the papers quoted in the table and 
assuming a solar Galactocentric distance of 8.5 kpc, consistent 
with Friel et al. 

In Col. 3 we list the ages of the sample clusters. Since they 
are retrieved from different sources, a unique and/or consistent 
age determination for them is not available. As is well known, 
age is very sensitive to the method employed, to the adopted 
stellar models in the case of isochrone fitting, and it can vary 
significantly from one study to another of the same cluster. In 
order to obtain a homogeneous set of ages for the oldest clus- 
ters (ages > 0.5 Gyr), we decided to use the morphological age 
indicator <SV (Phelps et al. 119941 ) and the metallicity-dependent 
calibration of Salaris et al. ( 120041 ). For most of the sample clus- 
ters £V is available from the literature; when not available, we 
obtained this quantity from available color-magnitude diagrams. 
Alternatively, we determined (5(B-V) and derived <SV from it 
(see again Phelps et al.). Note that for the calibrating clusters of 
Salaris et al. (|2p04) we directly used their age determination. For 
the younger clusters (the first seven lines of Table[TJ, we used the 
most recent age determinations available in the literature, such 
as those obtained with the lithium depletion boundary method. 
Use of non-homogeneous age determinations for young clusters 
will not affect our results and discussion, given our classification 
of the OCs in three major age bins (see below). 

Information on the observations and the analysis techniques 
is listed in Cols. 5, 6, 7, 8: specifically, we give the signal 
to noise ratio (S/N), the resolving power (R), the number and 
type of stars analyzed (giant or dwarf stars), and the method 
of analysis. The values of [Fe/H] and their errors, references 
to the observations and abundance analysis are listed in Cols. 9 
and 10. As the table shows, the values of R and S/N cover a 
rather wide interval, although for most samples the resolution 
is R~30,000-40,000. The metallicities are more often derived 
through equivalent width measurements, although in a few cases 
spectral synthesis was employed. The analysis is based on gi- 
ant stars for most of the intermediate age and old clusters, while 
for the young, more close-by cluster main sequence stars have 
been used. As discussed by Randich et al. (2006) for M67 and 
Pasquini et al. (2004) for IC4651, this should not introduce any 
major bias for the elements considered here. The numbers of 
stars used in each analysis range from 1 to 15, with the excep- 
tion of the very rich set by Paulson et al. (2003) for the Hyades 
(55 dwarf stars). 

In order to check whether (and to what extent) use of differ- 
ent abundance scales might introduce spurious trends or increase 
the dispersion, we compared different abundance measurements 
for the same cluster. Considering only the most recent studies 
based on high resolution, we found 11 cases with more than 
one [Fe/H] determination available. The difference between the 
[Fe/H] values does not depend on [Fe/H] itself and the average 
is rather small; namely, A[Fe/H]=0.12 dex with a standard de- 
viation of 0.08 dex. This value is somewhat larger (a factor of 
~ 2) than the typical error in [Fe/H] of most clusters; however, 
we believe that it is small enough not to affect our results in a 
significant way. For other element abundances, very few clusters 
have more than one determination; these studies indicate that 
the differences in [X/H] are likely larger and can be as large as 
0.2 dex and might introduce a certain amount of scatter in the 
observed distributions. 



3. The chemical evolution model 

The model adopted in this work is a generalization of the multi- 
phase model by Ferrini et al. ( 1992), originally built for the solar 
neighborhood, and subsequently extended to the entire Galaxy 
(Ferrini et al. [1994), and to other disk galaxies (e.g., Molla et 
al. [19961 [19971 Molla & Diaz [20051 Magrini et al. 120071 We 
refer to Ferrini et al ( 11992119941 ), Molla et al. ( fl996l ), and to 
Magrini et al. (2007) for a detailed description of the general 
model. The main assumptions of our model are that the Galaxy 
disk is formed by infall of gas from the halo and from the inter- 
stellar medium. The adopted infall has an exponentially decreas- 
ing law. This produces an inside-out formation scenario where 
the inner parts of the disk evolve more rapidly than the outer 
ones. We consider also an alternative parametrization of the in- 
fall including a constant amount of gas per unit area in addition 
to the exponentially decreasing law. This latter parametrization 
is useful in reproducing the metallicity gradient, but resulted in 
an enhanced star formation rate in the outer regions. 

In the Appendix the overall content of our model and in par- 
ticular the specific assumptions and parameterizations we used 
to reproduce the Galactic disk properties are described in detail. 

4. Radial gradients: comparison between models 
and observations 

In Figs. [1] [3] |U and |5j we compare abundances in OCs with 
the predictions of the model described in the Appendix for 
the a-elements Si, Ca, Ti, and the iron-peak elements, Fe, Ni, 
and Q0. In general, the observed abundance distributions show 
a two-slope feature with a change from a steep gradient for 
^gc ^ 11-12 kpc and a flat plateau at greater distances, at 
least in the two older age bins. For simplicity, we will assume 
12 kpc as the threshold Galactocentric distance to distinguish the 
inner regions from the outer plateau. We caution that this value is 
rather arbitrary since the samples in each age bin are limited and 
there is a lack of OC data in the radial range 12-15 kpc which 
would allow us to better constrain this value. We note however 
that our conclusions would not change by assuming a smaller 
(e.g., 1 1 kpc) or larger (e.g., 13 kpc) radius. For the Cepheids a 
change of slope also is seen at 7 kpc, with a steepening of the 
slope for Rqq < 7 kpc. 

In the following we will discuss separately the iron and other 
element gradients. For a better comparison with the chemical 
model, we have divided the cluster ages in three bins, namely (f) 
age less than or equal to 0.8 Gyr, (if) 0.8 Gyr<age< 4 Gyr, and 
(Hi) 4 Gyr<age< 1 1 Gyr. 

4.1. Iron 

In Fig. [T]and [2] we show the distribution of [Fe/H] of OCs and 
Cepheids as a function of the Galactocentric distance Rqc, and 
compare it with the [Fe/H] behaviour predicted by our model. 
For the Fe abundances of the Cepheids we have used the gradient 
derived by Andrievsky et al. (2004) and Lemasle et al. ( I2007I ). 
We recall that the Cepheids represent a young population; thus 
their data can be compared only with the current model gradient, 
while OC data allow a comparison with the model predictions at 
different epochs. 

The shape of the Fe gradient of OCs, at least in the most 
populated age bin (age between 4 and 1 1 Gyr) is characterized 

1 For the classification of elements in the a- and iron-peak we have 
followed Reddy et al. f2006b . 
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Fig. 1. Gradient of [Fe/H]: comparison between high-resolution 
data of OCs with the model decribed in the Appendix. Different 
panels show abundances for different age bins: top (panel 1), 
ages <0.8 Gyr; middle (panel 2), ages 0.8 Gyr < age < 4 Gyr; 
bottom (panel 3), ages 4 Gyr < age < 11 Gyr. The models 
adopted to compare with observations are for and 1 Gyr ago 
(panel 1), 1 and 4 Gyr ago (panel 2), and 4, 6, 8, 10, 11 Gyr ago 
(panel 3). 
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Fig. 2. Gradient of [Fe/H]: comparison between Cepheid data 
(Andrievsky et al. 2004) and our current model. 



by a negative slope at short radii and a change of slope at around 
~ 1 1 - 12 kpc, resulting in a plateau at larger distances. This be- 
haviour is qualitatively similar to that of Cepheid data, although 
the latter cover an inner part of the disk and show a change of 
slope around 7 kpc, but do not extend beyond 15 kpc. 

More specifically, it is useful to analyze in detail the gradient 
in the three radial zones selected, namely 4 kpc < Rcc S= 7 kpc 
(region I), 7 kpc < Rcc 55 12 kpc (region II), 12 Kpc < Rcc < 
22 kpc (region III). From Tableland Fig.|2]it is evident that the 
current gradient outlined by the metallicity of Cepheids is very 
steep in region I, where, unfortunately, no high-resolution OC 
data are available. In region II, the gradient traced by Cepheids 
(-0.051 and -0.044 dex kpc -1 according to Lemasle et al. 2007 
and Andrievsky et al. 120041 respectively) is consistent with that 
outlined by the youngest OCs (-0.053 dex kpc 1 -see Sect. [6j. 
Finally, as already suggested by the most recent papers (e.g. 
Yong et al. 120531 Carraro et al. |2004ll2007l Sestito et al.|2008), 
data from of OCs in region III indicate a very flat outer gradient, 
with a slope consistent with zero at all epochs (see also Sect.[6). 

As noted by Andrievsky et al. (2004), the existence of a dis- 
continuity in the metallicity gradient at ~ 11-12 kpc, observed 
for both Cepheid and OCs data, is predicted and explained by 
two chemical evolution models. One of them is the model pro- 
posed by Andrievsky et al. (2004) themselves, and the other is 
the model C presented by Chiappini et al. (2001). In particular, 
Andrievsky et al. (2004) assume that the SFR is a combination of 
two components: one proportional to the gas surface density, and 
the other depending upon the relative velocity of the interstellar 
gas and spiral arms. However, their model predicts a strong de- 
crease in metallicity beyond ~ 14 kpc, whereas the observed 
gradient from OCs is flat up to ~ 22 kpc. On the other hand, 
Chiappini et al. (2001) assume two main accretion episodes in 
the lifetime of the Galaxy, the first one forming the halo and 
bulge and the second one forming the thin disk. Both episodes 
are described by exponential laws. Their model C does not as- 
sume any threshold in the gas density during the halo/thick-disk 
phase, and thus allows the formation of the outer plateau from 
infalling gas enriched in the halo. 

Figs. Q]and[2] show that our model is also able to reproduce 
the flat gradient at large Galactocentric distances. This is ex- 
plained in the framework of our chemical evolution model by 
the radial dependence of the infall rate (exponentially decreas- 
ing with radius) which, combined with the gas distribution in the 
halo, results in a gas falling into the disk rapidly in the regions 
with R<12 kpc and at a almost constant rate in the outer regions 
(see Fig. IA.3I ). In addition, the radial dependence of the star and 
cloud formation processes, that are more efficient in the inner 
regions, produces a higher SFR in the inner regions, resulting 
in a higher metal production. In contrast, in the outer disk, the 
low metallicity of the infalling gas and the low SFR contribute 
to maintain a flat and slowly evolving metallicity gradient. 

To reproduce a completely flat gradient in the external re- 
gions, as indicated by the older OCs, our model would need 
an additional accretion of material uniformly falling onto the 
disk. This however would result in an inconsistent behaviour of 
the current SFR, predicting higher SFR in the external regions. 
We then considered several possible events that could explain 
the trend of the metallicity gradient traced by older OCs. The 
plateau is most evident only in the older group of OCs, and thus 
it could be related to peculiar episodes in the past history of the 
Galaxy. The outer plateau could be a signature of a past merger 
at high radii with pre-enriched material, or similarly to what is 
suggested by the model of Chiappini et al. (2001), of a gas halo 
already enriched. In this case it would not be necessary to mod- 
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ify the radial profile of the SFR to have higher element abun- 
dance at high radii, since the inflow of material would bring a 
higher content in metals. Alternatively, it could be caused by a 
strong episodic merger in the past with clouds/galaxies at high 
Galactocentric radii. Sporadic episodes of mergers with galaxies 
at large Galactocentric radii, more frequent in the past history of 
our Galaxy, would not strongly affect the current SFR and would 
explain the outer plateau in metallicity. However, fine tuning of 
the composition of the merger would be needed to reproduce the 
exact plateau value in iron and the distribution of other elements 
(see next section). 



4.2. Iron-peak- and a-elements 

Two iron-peak elements, Ni and Cr, and three a-elements (Si, 
Ca, and Ti) have been measured in a large number of OCs be- 
longing to the selected sample. Their gradients are plotted in 
Figs. [3] |4] and [5] The radial gradients are in most cases very 
similar and also very close to the iron gradient, with a nega- 
tive slope at low radii and a slightly decreasing distribution at 
larger distances (see Table |3). Note that in a few cases we were 
not able to derive a quantitative value for the slope, due to low 
statistics and/or to the small range of Rgc covered by the clus- 
ters with available abundances. Also, in two of cases (Ca and 
Cr) the slope in the 7-12 kpc range in the intermediate-age bin 
is consistent with 0, at variance with that of Fe. For all elements, 
however, the slope in the oldest age bin in the 7-12 kpc interval 
is very similar to the slope of [Fe/H] . This similarity is expected 
in our model because the metallicity gradient is determined for 
all elements by the radial dependence of the SFR and the infall. 

In order to reproduce the absolute mean values of 12 + 
log(Cr/H), we have assumed a production of Cr from SN la 
lower by a factor of ~ 2 than that by Iwamoto et al. (1999), 
while [Ni/H], [Ca/H], [Si/H], and [Ti/H] values are all we ll re- 
produced using the original yields by Iwamoto et al. (1999 ). The 
best agreement of the model prediction with the observed time 
and radial behaviours of the Ti abundance is obtained adopt- 
ing stellar yields from massive stars higher by a factor of ~ 2 
than those computed by Chieffi & Limongi (2004). This correc- 
tion is in agreement with the empirical stellar yields adopted by 
Francois et al. (|2004| l. 

In summary, all elements under consideration, in particular 
in the most represented age range, 4 Gyr < age < 10 Gyr, con- 
firm the existence of a steep gradient in the inner regions and a 
change of slope at around ~12 kpc with a flattening of the gradi- 
ent towards the outer regions. 



5. [X/Fe] radial gradients 

Although abundance ratios are less model-dependent than abso- 
lute abundances (see e.g. Matteucci & Chiappini 2005), never- 
theless they are affected by assumptions on stellar nucleosynthe- 
sis processes, the IMF, and stellar lifetimes. In particular, abun- 
dance ratios involving elements produced on short timescales by 
massive stars dying as SN II, like <?-elements, and elements pro- 
duced on longer timescales by e.g. SN I, like iron-peak elements, 
allow us to study the role of different types of supernovae in the 
chemical history of the Galaxy. 

In Fig. [6] [7] and|8] we compare the values [X/Fe] from high- 
resolution data with the predictions of the model. The production 
of iron-peak elements is dominated by SN la. Since iron has the 
same SN la origin, the ratios [X/Fe] of these elements are not 
expected to show any evident temporal evolution. 
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Fig. 3. Radial gradients of 12 + log(Cr/H) and 12 + log(Ni/Fe): 
comparison between high-resolution data of OCs with our 
model. Symbols and line-types are as in Fig. Q] These figures 
are available in electronic form only. 
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Fig. 5. Radial gradient of 12 + log(Ti/H): comparison between 
high-resolution data of OCs with our model decribed. Symbols 
and line-types are as in Fig. Q] These figures are available in 
electronic form only. 



Two iron-peak elements, Ni and Cr, have been measured in a 
large number of clusters belonging to our sample. Both elements 
show an almost constant average value, close to the solar one at 
all ages. The absolute mean values of [Cr/Fe] are reproduced 
assuming a production of Cr from SN la lower by a factor of ~ 2 
than that by Iwamoto et al. ( 1999 ), while [Ni/Fe] values are well 
reproduced using the original yields by Iwamoto et al. d 1 999b . 

Models and observations of both elements, in particular in 
the most represented age range, 4 Gyr < age < 1 1 Gyr, con- 
firm a constant behaviour of [X/Fe] with Galactocentric radius. 
We recall that typical formal uncertainties in [X/Fe] vary from 
0.01 to 0.15 dex. In the case of Ni lower values of [Ni/Fe] are 
observed in the few young clusters analyzed than in the older 
cluster sample. Also the [Cr/Fe] of some clusters in the old-age 
bin have lower values than expected from the model. Both these 
behaviours are not reproduced by the model, and they could be 
statistical variations, not significant, due to the small sample of 
clusters with these measurements. 

Data for the three c-elements observed in the high-resolution 
sample of OCs (Si, Ca, and Ti) indicate a rather flat distribu- 
tions with Galactocentric distance, and nearly solar [a/Fe] val- 
ues, slightly higher for older clusters. The largest contribution to 
the production of Si, Ca, and Ti comes from massive stars with 
M > 8M Q (e.g., Chieffi & Limongi |2004| ) while small amounts 
of these elements are also produced by SN la (e.g., Iwamoto et 
al. |19991 >. This is especially true for Ca: for example, SN la mod- 
els by Iwamoto et al. 0999) predict a typical synthesized mass 
of Ca of about 0.02-0.03 M Q per SN la, to be compared with the 
0.016 M e of Ca produced by a star of 25 M of solar metallic- 
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Table 2. Average [Fe/H] from OC data of different epochs and 
radial regions (see the text for the definition of radial regions). 



age (Gyr) 



Region II 



Region III 



< 0.8 
0.8-4 
4-10 



0.016 ±0.14 
0.05 ±0.11 
-0.04 ± 0.29 



-0.17 ±0.12 
-0.29 ±0.16 
-0.33 ±0.13 



ity (cf. Chieffi & Limongi 2004). Thus the production of Ca and 
Fe, integrated over the IMF, occur on comparable time scales, 
and the [Ca/Fe] gradient does not show any evident temporal 
evolution. Silicon, on the other hand, being mostly produced by 
rapidly evolving massive stars, shows high values compared to 
Fe in the past, and a more rapidly increasing [Si/Fe] radial gradi- 
ent than e.g., the [Ca/Fe] gradient, as suggested by e.g., Cescutti 
et al. (120061) . These temporal and spatial behaviours are well re- 
produced by the model, that also predicts a slight enhancement 
of [a/Fe] in the outer regions, due to the different spatial loca- 
tion of SN II and SN la. This is partially confirmed by the limited 
number of high-resolution observations of OCs. 

Titanium shows a large dispersion with cluster age, similarly 
to the one observed for [Si/Fe]. The best agreement of the model 
prediction with the observed temporal and radial behaviours of 
the Ti abundance is obtained adopting stellar yields from mas- 
sive stars higher by a factor of ~ 2 than those computed by 
Chieffi & Limongi (2004). This correction is in agreement with 
the empirical stellar yields adopted by Francois et al. (2004). 

As described in the Appendix, with our adopted model pa- 
rameters (infall, star and cloud formation efficiencies, chemi- 
cal yields, etc.) the best agreement of the theoretical gradients 
with the observations is obtained using the IMF by Kroupa et 
al. 0993). We also checked several other parameterizations of 
the IMF: assuming the "best" values (i.e. measured indepen- 
dently) for the various efficiencies of star and cloud formation 
at solar radii, we have varied the adopted parametrization of the 
IMF. A good agreement was also obtained with the IMFs by 
Ferrini et al. (|T990) and Scalo (1986). These parameterizations 
are characterized by a lower number of massive stars when com- 
pared e.g., w ith the Salpeter ( 119551 ), Scalo ( 119981 ), and Chabrier 
et al. (120031 ) IMFs. The high-resolution OC abundances of a- 
elements seem to weaken the use of these other parameteriza- 
tions that in fact result in an overproduction of massive stars, 
and, consequently, a production of a-elements greater than ob- 
served. For example, a Salpeter IMF would result in an over- 
production of Ti in the model of ~ 0.3 dex with respect to the 
observations. 



6. The time evolution of the metallicity gradient 

In their paper, Maciel et al. (2006) showed the results on the 
time evolution of the [Fe/H] gradient obtained by several au- 
thors, both from models (Hou et al. 120001 Chiappini et al. 120011 ) 
and from observations (PNe by Maciel et al. 2006, OCs by Friel 
et al. 2002, Chen et al. 2003, and H n regions, Cepheids, OB 
stars from several authors, see references therein). The metal- 
licity gradient was computed over the whole disk. However, as 
already remarked, the gradient does not have a constant slope 
over the whole disk, and, thus, a direct comparison with our data 
is difficult to obtain. 

Therefore we compared our results with the samples of OCs 
by Friel et al. (2002) and Chen et al. (2003) in the same radial 
regions. Their samples were discussed in the previous sections. 
In summary, they differ from our sample in data quality (high- 
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Fig. 6. Radial gradients of [Ni/Fe] and [Cr/Fe]: comparison be- 
tween high-resolution data of OCs with our model. Different 
panels show abundances for different age bins: top (panel 1), 
ages <0.8 Gyr; middle (panel 2), ages 0.8 Gyr < age < 4 Gyr; 
bottom (panel 3), ages 4 Gyr < age < 11 Gyr. The models 
adopted to compare with observations are for and 1 Gyr ago 
(panel 1), 1 and 4 Gyr ago (panel 2), and 4, 11 Gyr ago (panel 
3). The shaded regions indicate the area between the two mod- 
els shown in each panel. These figures are available in electronic 
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tween high-resolution data of OCs with our model. Symbols and 
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resolution spectroscopy vs. low-resolution spectroscopy and/or 
photometry) and radial range (they extend up to ~16 kpc, while 
our sample extends up to 22 kpc). The OCs located at high 
Galactocentric radius are fundamental in determining the slope 
in the older and intermediate age bin. In order to compare with 
the results of Friel et al. (2002) and Chen et al. (2003) we have 
re-computed the slope of the older and intermediate age bins 
for Rcc ^ 16, obtaining a slope for the cumulative sample of 
-0.067 dex kpc 1 , and for only the old bin of -0.079 dex kpc -1 , 
in closer agreement with their results: -0.075 dex kpc -1 for the 
sample of Chen et al. (2003) with ages >0.8 Gyr and -0.075 dex 
kpc -1 for the sample of Friel et al. (2002) with ages >4 Gyr. 
This points out the importance of OCs at larger Galactocentric 
distances for a correct determination of the outer gradient. 

The comparison with PNe is more tricky. The accuracy of 
the Galactic metallicity gradient derived for PNe has been de- 
bated for a long time because of the large uncertainties on PN 
distances. Recent results (Stanghellini et al. 120061 Perinotto & 
Morbidelli 120061 1 have shown very different gradients for the 
Galactic PN population than previous results. At variance with 
Maciel et al.'s results, Stanghellini et al. (2006) and Perinotto & 
Morbidelli (2006) found an almost flat gradient at all epochs. On 
the other hand, some other recent works find steeper gradients 
(e.g., Pottasch & Bernard-Salas 2006. Maciel & Costa [2008b 
pointing to the need for an accurate determination of the age in 
order to measure gradients. The main differences among these 
works are the adopted distance scales, demonstrating how this 
choice can affect the estimate of the metallicity gradient, in con- 
junction with different age estimates for the PN progenitors. 



Magrini et al,: Evolution of the Galactic metallicity gradient 



9 



Thus a faithful determination of the Galactic metallicity gradient 
from PNe is far from being obtained. 

A more appropriate way to compare the chemical evolution 
model with observations is to divide the Galaxy in several radial 
regions. In Tables [2] and [7] we list the the average abundances 
of iron in each age bin and radial region, and the best-fit slopes 
in the three radial regions. The latter were computed through a 
weighted fit procedure for the OCs and taken from Andrievsky 
et al. (120041 and Lemasle et al. (|2007T > for the Cepheids. The 
slopes of the models are also shown. As anticipated above, three 
age bins were considered: age less than or equal to 0.8 Gyr, in 
the first panel; 0.8 < age < 4 Gyr, in the second panel; and 
4 < age < 1 1 Gyr in the third panel. For OCs with ages between 
4 Gyr and 1 1 Gyr, we excluded from the calculation of gradient 
the metal-rich cluster NGC 679 1 (marked with an empty circle 
in Figs. Q] [3] I4I5K because it was probably born in a different 
place in the disk or could even have been captured by our Galaxy 
during a merger event (Carraro et al. 2006 but see also Bedin 
et al. 120061 for a contrasting view). Thus the age range for our 
further analysis is between 4 Gyr and 10 Gyr. 

First, we note that, although it is well known that a strict 
age-metallicity relationship does not hold in the Galactic disk, a 
tendency is present for older OC to have on average lower metal- 
licities. Second, each radial region shows a different behaviour 
due to the distinct ratios of SFR and infall. To better analyze the 
time evolution of the metallicity gradients, we plot in Fig.|9]the 
slopes of the gradient for the specific case of iron, the best stud- 
ied element, in the three radial regions defined above. Similar 
behaviours are observed also for the other elements analyzed in 
this study. 

As there are no observations of OCs in region I, it is not pos- 
sible to obtain the time evolution of the gradient in this region. 
In region II, our sample of OCs does not show any conclusive 
evidence for flattening with time of the slope of the gradient, al- 
though the present time gradient seems slightly flatter then the 
gradient in the two older age bins. Finally, in region III, the OC 
metallicities produce a gradient consistent with a zero slope at 
all epochs. 

Our data have allowed us to consider the time evolution in 
different regions of the disk, showing that in the ~ 7-22 kpc 
interval no marked evolution is present. 

The observations of Cepheids by Andrievsky et al. (2004) 
(see Tableland filled squares in Fig.|9]l show a good agreement 
with the absolute [Fe/H] predicted by the model in regions I, II, 
and III at the present time (see Table In region I, the model 
predicts a flattening of the metallicity gradient from the past 
epochs to the present time that cannot be tested because only 
Cepheid data are available in this region; in regions II and III 
it predicts an overall increase of metallicity at all radii but no 
significant change in the slope of the gradient (see Fig.|9]l. 

This can be understood as follows. In our model, the inner 
part of the disk is built by the rapid collapse of the halo, whereas 
the outer parts by the continuous slow accretion of intergalactic 
gas. In the inner disk, a steep chemical gradient is initially estab- 
lished by the strong concentration of infalling gas (see eq. IA.7K 
on a timescale probably of the order of a few Gyr or less (OCs 
are not able to constrain the chemical evolution of the gradient 
for the first ~ 5 Gyr of evolution of the Galaxy). The model 
predictions indicate that this initial gradient is preserved in the 
subsequent evolution of the disk. 

Turbulent mixing of the interstellar gas seems to be ineffi- 
cient in smoothing the initial metallicity gradient over the life- 
time of the Galaxy, in agreement with theoretical estimates; for 
example, according to Talbot & Arnett ( 119731 1, the timescale for 



mixing a region of a radial extent of 5 kpc like Region II is about 
75 Gyr. Also the role of large-scale radial flows, often invoked 
in the chemical evolution of disk galaxies, seems quite marginal, 
since their effect is to steepen the metallicity gradient (Lacey & 
Fall [19851 Tosi [19881 Portinari & Chio si EUOOl contrary to ob- 
servational evidence. However Friedli (1 19981 1 found that flows 
induced by bars can flatten the gradient. In the outer disk, on the 
other hand, accretion of intergalactic gas and disk evolution oc- 
cur almost independently of radius, and the metallicity gradient 
remains flat. 

We note that the slope of older OCs is significantly affected 
by the presence of a number of old clusters located at low 
Galactocentric distances with solar or super-solar metallicities. 
Old clusters have had enough time to move along their orbits, 
and the possibility exists that they were actually born in inner 
regions far from their present location. Although it seems un- 
likely that all old clusters with (over)-solar metallicities have 
different birthplaces, we stress that knowledge of the orbits of 
these older OCs, and thus their previous location, would help 
in definitively settling the issue. Also, a few clusters located at 
distances of 11-12 kpc are present with [Fe/H] somewhat be- 
low the plateau value observed at larger Galactocentric distances 
and characterized by large error bars. New observations for these 
clusters should be obtained in order to better constrain the slope 
of the gradient. Finally, metallicity determinations should be per- 
formed for both young and old clusters located at Galactocentric 
distances lower than ~ 7 kpc. 

7. Summary 

s We have assembled a homogeneous sample of abundances in 
OCs from high-resolution spectroscopic data and compared it to 
the results of a GCE model. We examined in particular the radial 
gradients of Fe/H, Si/H, Ca/H, Ti/H, Ni/H, and Cr/H, and also 
their [X/Fe] gradients at different epochs during the lifetime of 
the Galaxy. We have divided our data into three age bins and 
three radial regions. Our main results are the following: 

(i) From the observations of OCs and Cepheids (see Table [3), 
we find that the metallicity gradients of all the elements an- 
alyzed are very steep at low Galactocentric radii, in particu- 
lar for Rqc ^ 7 kpc. The gradients show a change of slope 
around 7 kpc, and further flattening at about 1 1-12 kpc, re- 
sulting in a metallicity plateau in the outer regions. 

(it) The combination of the observations of Cepheids and OCs 
allow us to follow the time evolution of the metallicity gra- 
dient and to compare it with the results of a GCE model. 
The time evolution of the gradient can be followed where 
observations of OCs with different ages are available for 
Rcc > 7 kpc. At Galactocentric radii 7 kpc < Rqc $12 kpc 
the combined observations of Cepheids and OCs suggest that 
the slope of the metallicity gradient has not changed sig- 
nificantly over the past ~ 10 Gyr; if anything, the gradient 
has flattened with time, as also suggested by the model for 
Rgc < 7 kpc. For 12 kpc < Rqc < 22 kpc, OCs belonging 
to different epochs and Cepheids all suggest a flat metallicity 
distribution. 

(Hi) The GCE model that better reproduces the data is character- 
ized by an inside-out formation of the disk. The infall of gas 
is represented by an exponential law which, combined with 
the distribution of gas in the halo, produces a rapid collapse 
in the regions with R<12 kpc and a uniform accretion in the 
outer regions. The inner regions are rapidly evolving due to 
the higher infall and SFR, while the outer parts evolve more 
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Fig. 9. Time evolution of the slope of the iron gradient in the 
three radial region analyzed: 4 kpc < Rgc ^ 7 kpc (region I), 
7 kpc < Rgc $12 kpc (region II), 12 kpc < Rgc ^ 22 kpc (re- 
gion III). Continuous lines are the slopes of the model gradients 
approximated with a linear fit. Filled squares are the slopes of 
the gradients of Cepheids. Filled circles represent the slopes of 
OC gradients with their uncertainties in the three age bins, as in 
Tabled 



slowly. The metallicity in the outer disk increases uniformly 
at all radii without a marked change of slope of the gradi- 
ent. Peculiar episodes of mergers in the past history of the 
Galaxy would explain the outer metallicity plateau indicated 
by the older clusters. 



Table 1. Open clusters included in our sample. We list in Col. 1 the cluster name, in Col. 2 the cluster age, in Cols. 3 and 4 the 
Galactocentric radius and its reference (an asterisk indicate that we computed Rgc assuming the cluster distance given in the quoted 
reference). In Cols. 5, 6, 7, 8, 9, and 10 we list the main properties of the observations and of the analysis: the signal to noise ratio, 
the resolving power, the number and type of stars analyzed, the adopted method (Equivalent widths, EW, or spectral synthesis), the 
iron abundance with its error, and the reference. 



Cluster 


Age 


^GC 


Ref. 


S/N 


R 


N.Stars 


Method 


[Fe/H] 


Ref. 




(Gyr) 


(kpc) 


(«GC) 












([Fe/H]) 



IC4665 

IC2602 

IC2391 

Blanco 1 

Pleiades 

NGC6475 

Mil 

M34 

NGC3960 

Hyades 

Praesepe 

NGC2324 

NGC2660 

Mel71 

Ru4 

Ru7 

NGC2360 

NGC2112 

NGC2477 

NGC7789 

NGC3680 

NGC6134 

IC4651 

Be73 

NGC2506 

NGC2141 

NGC6819 

Be66 

Be22 

M67 

Be20 

Be29 

NGC7142 

NGC6253 

NGC2243 

Be75 

Be31 

Mel66 

Be32 

Be25 

NGC188 

Saurerl 

Cr261 

Bel7 

NGC6791 



0.025 

0.035 

0.053 

0.1-0.15 

0.125 

0.220 

0.250 

0.250 

0.7 

0.7 

0.7 

0.7 

0.8 

0.8 

0.8 

0.8 

0.9 

1.0 

1.0 

1.1 

1.4 

1.6 

1.7 

1.9 

2.2 

2.5 

2.9 

3.7 

4.2 

4.3 

4.3 

4.3 

4.4 

4.5 

4.7 

5.1 

5.2 

5.4 

6.1 

6.1 

6.3 

6.6 

8.4 

10.07 

10.9 



8.18 

8.49 

8.50 

8.50 

8.60 

8.22 

6.86 

8.90 

7.96 

8.50 

8.62 

11.65 

9.18 

10.0 

12.40 

13.50 

9.28 

9.21 

8.94 

9.44 

8.27 

7.68 

7.70 

17.40 

10.88 

12.42 

8.18 

12.9 

11.92 

9.05 

16.37 

22.0 

9.70 

7.10 

10.76 

15.50 

12.02 

10.21 

11.35 

18.20 

9.35 

19.3 

7.52 

11.2 

8.12 



Manzi et al. 120081 * 

WEBDA*^ 

WEBDA - 

Friel (20061 

Friel 120061 

Robichon et al. fl999T 

Gonz alez & Wallerstein 120001 ' 

Friel (2006) 

Friel et al. (20021 

Friel 12006) 

Friel 179951 

Friel et al. 12002) 

Friel 1T995) 

Friel {2006) 

Carraro et al. 1 2007 1* 

Carraro et al. i 2007 l* 

Friel etal. 12002) 

Friel et al. {2002] 

Friel et al. (20021 

Friel et al. 12002) 

Friel et al. 12002) 

WEB DA* 

Friel 12006) 

Carraro etal. 1 2007 1 1 

Friel et al. 12TO) 

Friel et al. 12002) 

Friel et al. 12002) 

Phelp s & Ja nes fl996l ) 

Friel 11995) 

Friel et al. 12002) 

Friel etal.. 12002) 

Friel 12006) 

Jacobson et al. 12008) 

Bragaglia & Tosi (2006) 

Friel et al. 12002) 

Carraro et al. i 2007 l* 

Friel et al. 12002) 

Friel et al. {2002] 

Friel et al. (2002| 

Carraro et al. 1 2007 l* 

Friel et al. 12002) 

Friel 12006) 

Friel eTaT fe002] 

Friel et al. (20051 

Friel et al. 12002) 



30-150 


60,000 


18,dwa 


Synt 


-0.03 ± 0.04 


30-100 


18,000-44,000 


9,dwa 


EW 


-0.05 ± 0.05 


30-100 


18,000-44,000 


4,dwarfs 


EW 


-0.03 ± 0.07 


70-100 


50,000 


8,dwa 


EW 


+0.04 ± 0.04 


70 


45,000 


15,dwa 


EW 


-0.03 ± 0.02 


50-150 


29,000 


13,dwa 


EW 


+0.14 ±0.06 


85,130 


38,000 


10,gia 


EW/synthesis 


+0.17 ±0.13 


70 


45,000 


4,dwa 


EW 


+0.07 ± 0.03 


100-150 


45,000 


6,giants 


EW 


+0.02 ± 0.04 


100-200 


60,000 


55,dwa 


EW 


+0.13 ±0.05 


130 


100000 


7 


EW 


+0.27 ± 0.05 


90-150 


45,000 


7,gia 


EW 


-0.17 ±0.05 


50-80 


45,000 


5,gia 


EW 


+0.04 ± 0.04 


100 


16,000,34,000 


2,gia 


EW 


-0.32 ±0.16 


25^10 


40,000 


3,gia 


EW 


-0.09 ± 0.05 


25^10 


40,000 


5,gia 


EW 


-0.26 ± 0.05 


50-100 


28,000 


4,gia 


EW 


+0.07 ± 0.07 


70-80 


16,000,34,000 


2,gia 


EW 


+0.16 ±0.03 


100-150 


45,000 


6,gia 


EW 


+0.07 ± 0.03 


>50 


30,000 


9,gia 


EW/synthesis 


-0.04 ± 0.05 


80 


100,000 


2,dwa 


EW 


-0.04 ± 0.03 


80-200 


43,000 


3,gia 


EW 


+0.15 ±0.07 


100 


48,000 


4,gia 


EW 


+0.11 ±0.01 


25^10 


40,000 


2,gia 


EW 


-0.22 ±0.10 


80-110 


48,000 


3,gia 


EW 


-0.20 ± 0.02 


100 


28,000 


8,gia 


EW 


-0.18 ±0.15 


130 


40,000 


3,gia 


EW 


+0.09 ± 0.03 


5,15 


34,000 


2,gia 


EW 


-0.48 ± 0.24 


20-25 


34,000 


2,gia 


EW 


-0.32 ±0.19 


90-180 


45 000 


lU,UWd 


EW 


+0 0"* + 01 


40-80 


45^000 


2,gia 


EW 


-0.30 ± 0.02 


25-150 


45,000 


6,gia 


EW 


-0.31 ±0.03 


100 


30,000 


4,gia 


EW 


0.14 ±0.01 


60-140 


45,000 


5,gia 


EW 


+0.36 ± 0.07 


100 


30,000 


2,gia 


EW 


-0.48 ±0.15 


25^10 


40,000 


l.gia 


EW 


-0.22 ± 0.20 


100 


28,000 


5.gia 


EW 


-0.57 ± 0.23 


80-130 


45,000 


6,gia 


EW 


-0.33 ± 0.03 


50-100 


45,000 


9,gia 


EW 


-0.29 ± 0.04 


25-40 


40,000 


4,gia 


EW 


-0.20 ± 0.05 


20-35 


35,000-57,000 


5,dwa 


EW 


+0.01 ± 0.08 


80 


34,000 


2,gia 


EW 


-0.38 ±0.14 


70-130 


45,000 


7,gia 


EW 


+0.13 ±0.05 


15-100 


25,000 


4,gia 


EW 


-0.10 ±0.09 


10-60 


20,000 


15,gia 


synthesis 


+0.39 ± 0.08 



Shen etal. 120051 
Randich etal. (2001) 
Randich et al. (20 01) 
Ford et al. 12003) 
Boesgaard & Friel 11990) 
Sestito et al. 120031 
Gonzalez & Wallers tein (2000) 
Schuler et al. 12003) 
Sestito et al. 1 2006 i 
Paulson et al. i 2003 i 
Pace et al. 12008) 
Bragaglia et al. ( 2008 l 
Sestito et al. 120061 
Brown etal.(1 996) 
Carraro et al. i 2007 l 
Carraro et al. ( 2007 l 
Hamdani et al. i 2000 1 
Carrara et al. ( 2008 l 
Bragaglia et al. (20081 
Tautvaisiene et al. (20051 
Pace et al. 12008) 
Carretta et al. i 2004 1 
Carrettaet al.. ( 2004 1 
Carraro et al. i 2007 l 
Carretta et al. i 2004 1 
Yong et al. 120051 
Bragaglia et al. (2001) 
Villanova et al. {2005] 
Villanova et al . (20051 
Randich et al. ( 2006 1 
Sestito et al. 1 2008 i 
Sestito et al. 1 2008 i 
Jacobson et al. 1 2008 i 
Sestito et al. 120071 
Gratton & Co ntarini 119941 
Carraro et al. (20071 
Yong et al. (2005 i 
Sestito et al. 120081 
Sestito et al. 120061 
Carraro et al. ( 2007 l 
Randich et al. (2003) 
Carraro et al. ( 2004 1 
Sestito et al. 1 2008 i 
Friel et al. (20051 
Carraro et al. ( 2006 1 



2 

3 



< 
c 

o' 
3 



3" 

o 

a 



00 



Table 3. [Fe/H], [Ca/H], [Si/H], [Ti/H], [Ni/H], [Cr/H] gradients (dex kpcT 1 ) in different radial ranges. The observed gradients 
are computed from the data discussed in Sect. fusing a weighted linear fit. Model gradients are approximated by linear fits. T For 
OCs with 4 Gyr < age < 11 Gyr, the presence of the two metal rich OCs discussed in the text leads to a steep gradient (-0.16). 
Excluding them we still obtain flatter gradients (-0.13), but still steeper than model predictions. (fl) Andrievsky et al. (12)04) in the 
radial ranges: 4.0-6.6 kpc, 6.6-10.6 kpc, 10.6-14.6 kpc. w Lemasle et al. d2007l in the radial range 8-12 kpc. (e) OCs from Tableffl 
and references therein, (fl,) in the radial range 7-12 kpc 







model 






observations 








gradient 


range 


present time 


4 Gyr ago 


1 1 Gyr ago 


Cepheids 


OCs« r| 


OCs (fl 


OCs< c > 




(kpc) 








present time 


age < 0.8 Gyr 


0.8 Gyr < age < 4 Gyr 


4 Gyr < age < 1 1 Gyr 


A[Fe/H]/AS 


4-7 


-0.134 


-0.162 


-0.176 


-0.128 ±0.029'°' 


- 


- 


- 


A[Fe/H]/AS 


7-12 


-0.075 


-0.071 


-0.076 


-0.044 ± 0.01 l w 


-0.053 ± 0.029 


-0.094 ± 0.008 


-0.091 ±0.06 










-0.051±0.004 <D> 










A[Fe/H]/A« 


12-22 


-0.019 


-0.023 


-0.028 


+0.004 ± 0.01 1 (D » 


_ 


+0.005 ± 0.050 


-0.001 ±0.008 


A[Ca/H]/Afi 


4-7 


-0.137 


-0.172 


-0.156 




- 


- 


- 


A[Ca/H]/AR 


7-12 


-0.073 


-0.067 


-0.076 






0.02 ± 0.03 


-0.13 ±0.03 


A[Ca/H]/AR 


12-22 


-0.018 


-0.024 


-0.029 






-0.03 ± 0.09 


0.009 ± 0.009 


A[Si/H]/AR 


4-7 


-0.124 


-0.153 


-0.126 










A[Si/H]/AR 


7-12 


-0.063 


-0.056 


-0.076 




-0.05 ± 0.03 


-0.07 ±0.10 


-0.12 ±0.03 


A[Si/H]/A« 


12-22 


-0.017 


-0.024 


-0.029 








0.003 ±0.013 


A[Ti/H]/AS 


4-7 


-0.114 


-0.135 


-0.138 










A[Ti/H]/A« 


7-12 


-0.039 


-0.039 


-0.046 




-0.08 ± 0.03 


-0.07 ± 0.06 


-0.16 ±0.04 


A[Ti/H]/A« 


12-22 


-0.007 


-0.009 


-0.018 






-0.04 ± 0.02 


0.02 ± 0.01 


A[Ni/H]/AS 


4-7 


-0.126 


-0.157 


-0.156 










A[Ni/H]/A« 


7-12 


-0.072 


-0.071 


-0.076 






-0.07 ± 0.06 


-0.13 ±0.03 


A[Ni/H]/AS 


12-22 


-0.020 


-0.024 


-0.029 








0.00 ± 0.02 


A[Cr/H]/AS 


4-7 


-0.130 


-0.162 


-0.153 










A[Cr/H]/Afi 


7-12 


-0.072 


-0.068 


-0.075 




-0.08 ± 0.03 


0.04 ±0.11 


-0.18 ±0.03 


A[Cr/H]/AR 


12-22 


-0.020 


-0.024 


-0.029 








-0.006 ± 0.03 



Magrini et al.: Evolution of the Galactic metallicity gradient 



13 



Appendix A: The chemical evolution model of the 
Milky Way 

The model follows the evolution of the baryonic constituents 
(diffuse gas, clouds, stars, and remnants) of the Galaxy (halo 
and disk) from t — 0, when all the mass of the Galaxy is in the 
form of diffuse gas and it is located in the halo, until later times, 
when the mass fraction in the various components is modified 
by several conversion processes. 

The Galaxy is divided into N coaxial cylindrical annuli with 
inner and outer Galactocentric radii (i = l,N) and Ri+i, 
respectively, mean radius /?;+i/2 = (Ri + Ri+\)/2 and height 
h(Ri + \/2)' Each annulus is divided into two zones: the halo and 
the disk. These components are both made up of diffuse gas g, 
clouds c, stars s and stellar remnants r. The halo component is 
intended here quite generally as the material external to the disk 
from which it is formed and accreted, i.e. the primordial bary- 
onic halo, and/or the material accreted from interactions with 
small nearby dwarf galaxies or from the intergalactic medium 
during the lifetime of the Galaxy. 

The main processes considered in the model are: (J) conver- 
sion of diffuse gas into clouds; (z'z) formation of stars from clouds 
and (Hi) disruption of clouds by previous generations of massive 
stars; (iv) evolution of stars into remnants and return of a fraction 
of their mass to the diffuse gas. We describe in the following the 
main equations of the model, using subscripts H or D to indicate 
the halo or the disk, respectively. For a spherical halo of radius 
Rn+i = Rh, 



MRm/2) = (4 - Rm/t) V2 > 

and the halo volume in the i— th annulus is then 

V H ,i = it{RM-Rl)h(Ri+ii2). 



(A.l) 



(A.2) 



If Ro is the disk scale height (assumed independent of radius), 
the volume of the disk in each annulus is 



V DJ = n(R 2 M - Rj)R D . 



(A.3) 



At time t — the Galaxy is constituted by diffuse gas in the 
halo. At later times, the mass fraction in the various compo- 
nents is modified by several conversion processes: diffuse gas 
is converted into clouds, clouds collapse to form stars and are 
disrupted by massive stars, stars evolve into remnants and return 
a fraction of their mass to the diffuse gas. A disk of mass Mjyif) 
is formed by continuous infall from the halo of mass M#(f) at a 
rate 



dM L 
dt 



= fM H , 



(A.4) 



where / is a coefficient of the order of the inverse of the infall 
time scale. Clouds condense out of diffuse gas at a rate fj. and are 
disrupted by cloud-cloud collisions at a rate H', 



dM c ,,3/2 



H'Mi 



(A.5) 



where M g (t) and M c (t) are the mass fractions of diffuse gas and 
clouds, respectively. Stars form in the halo and the disk by cloud- 
cloud collisions at a rate H and by the interactions of massive 
stars with clouds at a rate a, 



dM s 
dt 



= HMt + aM s M c - DM S , 



(A.6) 



where M s (t) is the mass fraction in stars and D is the stellar death 
rate. 



Each annulus is then evolved independently (i.e. without ra- 
dial mass flows) keeping its total mass fixed from t = to 
fgai = 13.6 Gyr, the age of the Galaxy, computing the fraction 
of mass in each component in the two zones, and the chemical 
composition of the gas (assumed identical for diffuse gas and 
clouds). The rate coefficients of the model are all assumed to be 
independent of time but functions of the Galactocentric radius 
R. Their radial dependence is discussed in detail by Ferrini et 
al. ( fl994l l, Molla et al. ( fl996l and Molla & Diaz d20U5l In gen- 
eral, coefficients representing condensation processes (e.g. the 
formation of clouds from diffuse gas), being proportional to the 
inverse of the dynamical time, scale with the inverse square root 
of the zone volume, whereas the coefficients of binary processes 
(e.g. the formation of stars by cloud-cloud collisions) scale as 
the inverse of the zone volume; the coefficient of star formation 
induced by stars is independent of radius. 

A.1. The infall 

The accretion rate of gas in the disk E g , the infall law, is given in 
our model by 



3Lg(R) Z H (R) 

— ^ — = < 

dt t 



-t/T-R/A D 



(A.7) 



which represents the infall on the disk resulting from the col- 
lapse of the halo on a timescale t with radial lenghtscale Ap. For 
r and Ad we assume values typical for the solar neighborhood, 
t = 8 Gyr and Aq=2 kpc (Sackett 1997). The collapse timescale 
at the solar radius is based on several observational constraints 
(e.g., G-dwarf metallicity distribution, O/Fe vs. Fe/H relations, 
present-time infall rate), all suggesting a time scale for the for- 
mation of the disk at the solar circle of about 8 Gyr (see Molla 
& Diaz 2005 and references therein). 

We have also explored the consequence of an additional in- 
fall rate constant in time and uniform in radius. This is expressed 
by the following equation 



dZJR) 



^ V ' V = Zh(R) - t /T -RUn + I ^2 

dt t \ dt 



(A.8) 



where the first term is as in eq. (|A.7| > and the second term is 
represents continuous accretion of matter. This latter component, 
when tuned to the appropriate value of ~1 M Q yr~' integrated 
over 22 kpc radii, allows us to reproduce better the observed 
abundance plateau in the outer Galaxy up to Rqc ~ 22 kpc, as 
shown in the paper. On the other hand, the addition of a constant 
infall produces a radial profile of the SFR that is not consistent 
with current observations. In fact, in this case, we would have 
at higher Galactocentric distance, a SFR comparable to that at 
solar radius. For this reason, even if the outer flat gradient was 
well reproduced, we did not use this parametrization. 

A.2. The parameters of the model 

The mass fractions in each component computed by the model 
in each annulus are then converted into surface densities using as 
a normalization the observed total baryonic surface density (i.e. 
the sum of the Hi, H2, stellar and remnant surface densities). The 
observed gas surface density is taken from Dame et al. dl9931 l, 
and the stellar surface density, including remnants, from Boissier 
& Prantzos (119991 ), with an exponential scale-length of 2.5 kpc 
(Ruphy et al. 1996) and a normalization value at the solar radius 
of 43 + 5 M pc~ 2 (Mera et al. |T998l ). In the comparison with the 
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Fig.A.l. The radial profile of the SFR, normalized to its value 
at R Q =8 kpc and at the present time obtained with an exponen- 
tial infall law as in eq. ( 1A.71 i. The model curves are compared 
with the observations by Lacey & Fall (1985)— filled squares- 
, Gusten & Mezger (1983)— filled triangles- and Williams & 
McKee (1997)— filled circles. The model curves are at 1 (dotted- 
dashed line), 5 (dotted line), 8 (dashed line) Gyr from the forma- 
tion of the Galaxy, and the present time (continuous line). 



data, we identify the diffuse gas and cloud components with the 
Hi and the H2 gas, respectively. 

For our model of the Galaxy we adopt the "best" values 
of the various coefficients (infall, cloud and star formation ef- 
ficiencies) in the solar neighborhood, scaling them at different 
Galactocentric radii, as discussed by Ferrini et al. ( 1994), Molla 
et al. ([19961 Molla & Diaz d2005l Magrini et al. (120071 . The 
efficiencies of the processes regulating the conversion of diffuse 
gas into clouds and clouds into stars are close to those found in 
the best model for the Galaxy of Molla & Diaz (120051 1. corre- 
sponding to the efficiencies of their morphological type N - A 
(see their Table 2). In Table [ATI we show the values of the co- 
efficients we have adopted for our best model. Rq and Rh are 
the scale heights of the disk (cf. Sackett [T997l l and of the halo 
(cf. Boissier & Prantzos 1999), respectively. Ad is a typical disk 
scale length (cf. Ferrini et al. 1 19941 Portinari & Chiosi |1999T >. \i 
is the efficiency of the cloud formation from diffuse gas, H that 
of the star formation from cloud collisions, and H' that of cloud 
dispersion. 

The chemical enrichment of the gas is modeled using the 
matrix formalism developed by Talbot & Arnett dl9731 l. The el- 
ements of the restitution matrices Qjj(M, Z) are defined as the 
fraction of the mass of an element j initially present in a star 
of mass M and metallicity Z that it is converted into an el- 
ement i and ejected. Our GCE model takes into account two 
different metallicities, Z = 0.02 and Z = 0.006, and 22 stel- 
lar mass bins (21 for Z = 0.006) from M min = 0.8 M to 
M mdx = 100 M , for a total of 43 restitution matrices. For low- 
and intermediate-mass stars (M < 8 M ) we use the yields by 
Gavilan et al. (2005) for both values of the metallicity. For stars 
in the mass range 8 M < M < 35 M we adopt the yields by 
Chieffi & Limongi d2004l > for Z = 0.006 and Z = 0.02. 

We estimate the yields of stars in the mass range 35 M Q < 
M < 100 M B , which are not included in tables of Chieffi & 
Limongi (2004) by linear extrapolation of the yields in the mass 
range 8 M Q < M < 35 M . In the case of Ti, we increased the 
yields by a factor of 2 to reproduce the observations. The SNIa 
yields are taken from the model CDD1 by Iwamoto et al. ( 11999b . 



Another fundamental ingredient of a GCE model is the ini- 
tial mass function (IMF). Several parameterizations of the IMF 
have been widely used in the literature. Since the magnitude and 
the slope of the modeled chemical abundance gradients are re- 
lated to the choice of the IMF, we checked the consistency of 
several IMFs with the observed gradients. With the adopted in- 
fall, efficiencies and set of chemical yields, we have found the 
best agreement with the observed gradients using the IMF by 
Kroupa et al. ( [1993] ). 

A3. Results of the model 

A successful chemical evolution model of the Milky Way should 
reproduce the observed radial profiles and the integrated quan- 
tities. Fundamental radial profiles to be reproduced are: i) the 
gas radial distribution (e.g. Dame 1993); ii) the stellar exponen- 
tially decreasing profile; iii) the strongly decreasing SFR profile; 
iv) the metallicity gradients of various elements. The integrated 
quantities are instead the average infall rate, the average SFR, 
and the supernovae rates of both SNII and SNIa. 

In Figures IA.2I and IA.3I we present the comparison of the 
results of our model with the previously described radial profiles. 
Note that the model describes only the evolution of the disk, 
thus the inner regions hosting the bulge are not expected to be 
reproduced (R<5 kpc). The radial profiles of the total gas and of 
the stars at the present time are well reproduced, as is the radial 
profile of O/H as drawn by H 11 regions. Note however the large 
dispersion in the observations. The general trend of the SFR is 
reproduced by the model, even if the observed maximum is at 
~6-7 kpc, while the model maximum is located at ~9-10 kpc 

Our "best" model is also able to reproduce the observational 
integrated quantities for the Galaxy, as: (2) the present-time in- 
fall rate is 3.5 M yr~' integrated over 22 kpc radii, consistent 
with the recent estimate ~ 5 M Q yr _1 by de Boer (2004), and 
favoured by recent models of disk galaxy evolution (e.g., Naab 
& Ostriker 2006); (ii) the predicted current star formation rate 
is ~ 3 M Q yr integrated over 22 kpc radii, to be compared 
with the observed value e.g. by Pagel ( 1 19971 ) of 2-6 M yr -1 ; 
(iii) the current supernova rates obtained from the model are 
1.6 SN/century for SNIa and 1.4 SN/century for SNII, of the 
same order as the values for Sb galaxies ~ 1.8 SN/century and 
~ 2 SN/century, respectively (Mannucci et al. 2005). 
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